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Abstract 

The input to the configuration-interaction shell model includes many dozens or hundreds of inde- 
pendent two-body matrix elements. Previous studies have shown that when fitting to experimental 
low-lying spectra, the greatest sensitivity is to only a few linear combinations of matrix elements. 
Here we consider interactions drawn from the two-body random ensemble, or TBRE, and find that 
the low-lying spectra are also most sensitive to only a few linear combinations of two-body matrix 
elements, in a fashion nearly indistinguishable from an interaction empirically fit to data. We find 
in particular the spectra for both the random and empirical interactions are sensitive to similar 
matrix elements, which we analyze using monopole and contact interactions. 
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I. INTRODUCTION AND MOTIVATION 



The configuration-interaction shell model is a useful framework for a detailed under- 
standing of low-energy nuclear structure [l-3]. The many-body basis is a large dimension 
(10 3-10 ) set of Slater determinants, which are antisymmeterized products of single-particle 
states. One must truncate the single-particle states, corresponding to one or a few shells 
(typically using the harmonic oscillator as an approximation to the mean- field); the many- 
body basis may be further truncated. For phenomenological calculations one writes the 
Hamiltonian in terms of single-particle energies and two-body matrix elements, while for ab 
initio calculations one may extend this to three-body interactions {4]. 

The two-body matrix elements are the matrix elements of the residual interaction in the 
lab frame, 

V JT (ij,kl) = (ij;JT\V\kl;JT) (1) 

where \ij; JT) is the normalized, antisymmeterized product of particles in orbits labeled 
by i and j and coupled to good angular momentum J and isospin T. If one starts from a 
translationally invariant interaction between particles, one can either compute the integral 
in the lab frame or start in the relative frame and then transform to the lab frame; in either 
case there are correlations between the matrix elements, although they are not obvious to 
the casual observer. 

Often for semi-phenomenological calculations, one starts from a "realistic" interaction, 
and then adjusts the two-body matrix elements until the rms error on a set of experimentally 
known energy levels is minimized In the ls-Od shell, such a semi-phenomenological 

n n 

interaction has been recently derived pj, improving on an earlier interaction |6|]. 

It has been found that the fits are empirically dominated by a few linear combinations of 
matrix elements. The physical meaning of those dominant combinations is not immediately 
obvious. One might naively guess the linear correlations are due to an underlying transla- 
tionally invariant interaction (although a density dependence would destroy this). Somewhat 
more phenomenologically, it has been argued by appealing to mean-field properties that one 
can improve fits primarily through adjusting the monopole-monopole part of the interaction, 
that is, interaction terms that look like n a {ni > — 5 a b), where n a is the number of particles in 
the ath orbit. This protocol for shifting monopole strengths has been successfully applied 
to several semi-empirical interactions |7H11| 
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A related study 12|, investigating the origin of many-body forces from truncation of 
the model space, also found an empirical fit dominated by a few linear combinations of 
matrix elements. While much of the fit was dominated by the monopole interactions, even 
better agreement was brou ght about using a contact interaction motivated by its usage in 



mean-field calculations 
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14| and effective field theory 15j . 



In investigating the character and origin of the dominant matrix elements, it is useful to 
ask if there is anything special about the nuclear interaction. One way to ask this question 
is to compare with interactions drawn from the two-body random ensemble (TBRE), which 
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despite their arbitrary nature are known to echo some features of real nuclear spectra 

In this paper we conduct a sensitivity analysis of the low-lying spectra of random 
interactions and compare against a standard empirical interaction, USDB. We find that for 
more measures all the interactions are nearly indistinguishable, at least on a statistical level. 



II. METHODOLOGY AND RESULTS 

Our methodology follows previous work [l|, 0, [l2|; we work in the ls-Od valence space 
with an inert 16 core. Given an input set of two-body matrix elements (we leave aside 
single-particle energies and any A-dependent scaling), which we write as a vector v, we can 
calculate the eigenvalues E a {v) of the many-body Hamiltonian. For this work the label a 
ranged over all nuclides with < faience < A va i ence < 10 and took the ground state binding 
energy and the first five excitation energies. 

If one has a target spectrum E^, say from experiment, then the goal of the fit is to 
minimize 

E(^)-<) 2 (2) 

a 

(for simplicity we leave off the experimental uncertainty in each state). Expanding to first 
order 

dE 

E a (v + 5v) « E a (v) + J25v l ^- (3) 



then minimizing ([2]) yields 
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FIG. 1: (Color online)The spectra of eigenvalues from a singular-value decomposition of the sen- 
sitivity matrix A (Eq. [5]) for the two-body random ensemble. The lower curve is for ground state 
energies, while the upper curve is for excitation energies. 



This equation is in the form b = Ax where 



A, 



E 



dE a dE a 



a dvi dvj 



The derivatives come via the Hellmann-Feynman theorem 20 1 

dE a 



dvi 



H 



(5) 



(6) 



where Hi is the Hamitonian operator whose strength is t>j. 

We then find the eigenvalues of A (which is equivalent to finding the squares of the 
eigenvalues in the singular- value decomposition of dE a /dvi). We do this for both the USDB 
interaction and for an ensemble of 100 sets of random two-body interactions, also called the 
two-body random ensemble (TBRE). The results are shown for the TBRE in Fig. 1, where 
we have separated out the sensitivity just for the binding energies (ground state energies) 
and the excitation energies. Although not shown, the equivalent SVD eigenvalues for USDB 
are completely within the TBRE results. 

The lower curve is for ground states only, while the upper curve is for excitations energies 
relative to the ground state. Clearly, and perhaps unsurprisingly, the ground state energies 
are predominantly sensitive to just a few linear combinations of matrix elements-significantly 
fewer than excitations energies. 

Fig. 1 characterizes eigenvalues of A. The next step is to characterize the eigenvectors 
associated with the dominant eigenvalues, specifically by comparing with monopole and 
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FIG. 2: (Color online) The SVD spectrum that measures the overlap of the subspace defined by 
the six largest eigenvectors from Fig. 1, with the subspace defined by the monopole-monopole 
interaction. (Black) squares are for the TBRE, while (red) diamonds are for USDB. 

contact interactions. To do so we first discuss a method for quantifying the overlap between 
two vector subspaces 
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22J. Consider two vectors subspaces, Si and S , 2 . Let Vi be a 



matrix whose column vectors are the (orthonormal) basis vectors of Si, and similarly with 
V2. From these one constructs the overlap matrix ft = VjV2. Note that if the subspaces 
are not of equal dimension then O is not a square matrix. In any case we do a singular value 
decomposition of f2; the SVD eigenvalue spectrum then is a measure of the overlap of the 
two spaces. If the two spaces perfectly overlap then all eigenvalues are 1, if just N of the 
dimensions perfectly overlap than N eigenvalues will be 1 and the rest zero. Note that this 
method is invariant under arbitrary choice of orthonormal bases. 

We begin with the monopole-monopole interaction of the form n a (nb — S a b), which has 
six unique terms, and thus six vectors or linear combinations of matrix elements, in the 
s<i-shell. These we combine with the k most dominant linear combinations that arise from 
the previous analysis; somewhat arbitrarily we chose k = 6 (our results do not change 
qualitatively for other small values of k). The results, the SVD eigenvalues of Q, are shown 
in Fig. 2. It is important to note we compute the spectrum separately for each randomly 
generated interaction and then compute the distribution. 

The results for ground states and for excited states are similar, so we combine all states 
into a single calculation. The eigenvalues for USDB are roughly 50% higher than for the 
TBRE, but otherwise qualitatively very similar. 
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FIG. 3: SVD spectrum from the overlap of the dominant eigenvectors from USDB and the TBRE. 

We also compared for contact interactions; we took only two terms, the S = 0, T — 
1 channel and S — 1, T — channel (there being only the s-wave channel in relative 
coordinates). These results we summarize in Table 1. 



TABLE I: Leading eigenvalues from SVD of subspace overlaps of two-term contact interaction 



Interaction 


ground states excited states all states 


USDB 
TBRE 


0.60 0.58 0.62 
0.55 ±0.04 0.41 ±0.06 0.55 ± 0.04 



For comparison, the leading eigenvalue for the overlaps of USDB versus the six-term 
monopole is 0.94, while that of the TBRE versus monopole is 0.66 ± 0.03. There is some- 
what more sensitivity to the monopole interaction than the contact interaction; however, 
the reader should keep in mind that is not the whole story. Recall that when fitting an 
interaction, the linearized equations are cast in the form Ax = B. Our analysis in this 
paper is entirely with the eigenvalues of A, but in any fit one must also look to b (which in 
practice is the deviation of the theoretical spectra from experiment). For example, in 12] it 
was found that using a contact interaction brought better agreement than a monopole inter- 
action. One can understand this in terms of conjugate gradient methods 23J|: the direction 
of local steepest gradient may not in fact point towards the global minimum. 

By our measures so far, both the TBRE and the empirically-fit USDB look qualitatively 
similar. Therefore we take a final analysis by comparing the dominant linear combinations 
of the USDB with those from the TBRE. This is show in Fig. 3, using the same analysis as 
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for Fig. 2. For comparison with the previous results, the leading eigenvalue is 0.72 ± 0.06. 



III. CONCLUSION 

We have analyzed the sensitivity of the low-lying spectra of the random two-body ensem- 
ble of interactions to variations of the Hamiltonian matrix elements; by using singular value 
decomposition, we find the dominant linear combinations, which would be important in any 
fit to experimental data. We found the SVD eigenvalues follow a pattern remarkably simi- 
lar to that shown by semi- realistic/semi-phenomeno logical interactions such as USDB. We 
also analyzed the most dominant linear combinations of matrix elements by computing the 
overlap with monopole and contact interactions. Overall, both the TBRE and the empirical 
USDB had qualitatively similar results. 

The U.S. Department of Energy supported this investigation through contracts DE-FG02- 
96ER40985 and DE-FC02-09ER41587, and through subcontract B576152 by Lawrence Liv- 
ermore National Laboratory under contract DE-AC52-07NA27344. 
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